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Abstract 

We study the inclusive electron spectrum of ^ XuCi' decays near the kinematic limit for B — > 
Xcev transitions using a sample of 88 million T(4S') decays recorded by the BABAR detector at the 
PEP II e^e~- storage rings. For the electron momentum interval of 2.0 — 2.6GeV/c the partial 
branching fraction is measured to be AB{B X^eu) = (0.480 ± 0.029stat ± 0.053syst) x 10~^. 
Combining this result on A^B with measurements of the inclusive photon spectrum m. B ^ Xs'y 
transition we find \Vub\ = (3.94zb0.25expi0.37j^ zbO.lOj/ieorj;) x 10~^, where the first error represents 
the combined statistical and systematic experimental uncertainty of the partial branching fraction 
measurement, the second refers to the uncertainty of the determination of the fraction /„ of the 
inclusive electron spectrum that falls within our cuts, and the third error is due to theoretical 
uncertainties in the QCD corrections, plus the uncertainty in the 6-quark mass. All results are 
preliminary. 
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1 Introduction 



The increasingly precise measurements of CP asymmetries in i?-meson decays allow for stringent 
experimental tests of the Standard Model mechanism for CP violation [1] via the non-trivial phase 
in the Cabibbo-Kobayashi-Maskawa matrix. Improved measurements of CKM element iV^bl, the 
coupling of the b quark to the u quark, will enhance the sensitivity of such experimental tests. 

The extraction of \ Vub\ is a challenge, both theoretically and experimentally. Experimentally, the 
principal challenge is to separate the signal B X^eu decays from the 50 times larger B X^eu 
background. This is achieved by selecting regions of phase space in which the background is 
highly suppressed. In this paper we present a measurement of the inclusive electron spectrum for 
charmless semileptonic B decays, averaged over charged and neutral B mesons, near the kinematic 
endpoint. In the rest frame of the B meson, the kinematic endpoint of the electron spectrum for 
the dominant B —> XcCu decays is ~ 2.3GeV/c and ~ 2.6GeV/c for B —>■ X^ev decays. The 
finite momentum of the B mesons in the T{AS) decays causes additional spread of the electron 
momenta of ~ 200MeV/c, extending the endpoints to higher momenta. A narrow interval of 
about 300MeV/c remains dominated by electrons from B XuCU transitions, and this allows 
for a relatively precise measurement of the partial branching fractions in an interval that covers 
approximately 10-15% of the total electron spectrum for charmless semileptonic B decays free 
from significant BB background. In this analysis we extend the interval for signal extraction up to 
600 McV/ c covering ~ 30% of the total electron spectrum. 

Theoretically the weak decay rate for b —>■ ulu can be easily calculated at the parton level. It 
is proportional to \Vuh\^ and m^, where m\, refers to the b quark mass. To relate the B meson 
decay rate to \Vub\; the parton- level calculations have to be corrected for perturbative and non- 
perturbative QCD effects. These corrections can be calculated using operator-product expansions 
in powers of l/m;, and as [2]. However, near the kinematic endpoint in the lepton spectrum, 
these calculations break down because the spectra are affected by the "shape function", i.e., the 
distribution of the 6-quark momentum inside the B meson [3] , in addition to weak annihilation and 
other effects. Thus extrapolation from the limited momentum range near the endpoint to the full 
spectrum is a difficult task. At present, the shape function cannot be calculated, but it should be 
a universal property of the B meson. To leading order it must be the same for aWb ^ q transitions 
(here q represents any light quark). It has been proposed [4], [5] that \Vub\ can be extracted by 
comparing the high-energy end of the lepton spectrum with the high end of the photon spectrum 
in 6 — > S7 transitions, thus reducing the theoretical uncertainty on the shape function. Inclusive 
measurements of the photon spectrum are very challenging, but two results are now available [6] , [7] . 
This analysis is based on the same method as previous measurements of the lepton spectrum near 
the endpoint [8] (ARGUS), [9], [10](CLEO). The measurement of the partial branching fraction for 
charmless semileptonic decays and the extraction of \Vub\ presented here are updates of preliminary 
results [11] based on one-fourth the data. 

2 Data Sample, Detector, and Simulation 

The data used in this analysis were recorded with the BABMi detector at the PEP-II energy- 
asymmetric e+e~ collider. The data sample of 88.4 million BB events, corresponding to an in- 
tegrated luminosity of 80.4 fb~^, was collected at the T(4S') resonance; an additional sample of 

9.5 fb~^ was recorded at a center-of-mass energy just below the BB threshold. The second data 
set is used for subtraction of the ivoiy-BB contributions from the data collected on the T(4S') 
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resonance. The relative normalization of the two data samples has been derived from luminosity 
measurements, which are based on the number of detected /j.^/j.'^ pairs and the QED cross section 
for e+e~ production, adjusted for the small difference in center-of-mass energy. 

The BABAR detector has been described in detail elsewhere [12]. The most important com- 
ponents for this study are the charged-particle tracking system, consisting of a five-layer silicon 
detector and a 40-layer drift chamber, and the electromagnetic calorimeter assembled from 6580 
CsI(Tl) crystals. These detector components are embedded in a 1.5 T solenoidal magnetic field. 
Electron candidates are selected on the basis of the ratio of the energy detected in the calorimeter 
to the track momentum, the calorimeter shower shape, the energy loss in the drift chamber, and 
the angle reconstructed in the ring-imaging Cherenkov detector. 

The electron identification efficiency and the probabilities to misidentify a pion, kaon, or proton 
as an electron have been measured as a function of laboratory momentum and angles [13] with 
clean samples of tracks that were selected from data. This experimental information is used in 
the Monte Carlo simulation to improve the agreement with the data. Within the acceptance of 
the calorimeter, defined by the polar angle in the laboratory frame, —0.72 < cos 9iab < 0.92, the 
average electron efficiency is 92%. The average hadron misidentification rate is about 0.1%. 

Tracking efficiencies and resolution have been studied in great detail. Comparisons with the 
simulation have revealed small differences, which have been taken into account. No significant 
impact of non-Gaussian resolution tails has been found for high momentum tracks contributing to 
the endpoint region. 

We use Monte Carlo simulation of BB events to estimate signal efficiencies and background 
distributions. Most of the branching fractions for hadronic B and D decays and form factors are 
based on values reported in the Review of Particle Physics [14]. Charmless semileptonic decays, 
B — > Xuii/, are simulated as db mixture of thrcc-bociy decays (^X^ = TT,rj, p,ui, ...) based on the 
ISGW2 model [15], and decays to hadronic states Xu, with masses extending from 2mT^ to 
about 3.5 GeV/c^, according to the prescription of Ref. [16]. The hadronization of is performed 
by JETSET [17]. The motion of the b quark inside the B meson is implemented with the shape 
function parameterization given in [16]. The low-mass resonant decays are mixed with the non- 
resonant states in such a way as to retain the cumulative distribution in the hadron mass mxy, as 
predicted by the non-resonant model of [16]. 

For the simulation of the dominant B — > Xciv decays we have chosen a variety of models. We 
use ISGW2 [15] for B Diu and the various decays to higher mass D** resonances, adopt a 
prescription by Goity- Roberts [18] for the non-resonant B D^*^-k(!.v decays, and use an HQET 
parameterization [19] of the form factors for B D*ii'. For these form factors we use the recent 
results from BABMl [20] , which have much smaller uncertainties than the earlier measurements by 
CLEO [21]. The branching fractions for individual B Xciiy decays are adjusted to match the 
data (see below). 

The Monte Carlo simulations include radiative effects such as bremsstrahlung in the detector 
material and QED initial and final state radiation. Adjustments for small variation of the beam 
energy over time have also been included. 
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3 Analysis 



3.1 Event Selection 

We select semileptonic -B-decay events by requiring that there be an electron with momentum 
Pf, > 1.1 GeV/c in the T(45) rest frame. Throughout this paper, all quantities are measured in the 
T{AS) rest frame unless it is specified otherwise. To reject electrons from the decay J/V' e+e" we 
combine the electron candidate with any second electron of opposite charge and reject the candidate 
if the invariant mass of the pair falls in the interval 3.00 < ruee < 3.15GeV/c^. 

To suppress background from non-BB events, primarily low multiplicity QED (including r+r" 
pairs) and e"'"e~ — > qq processes, we veto events with less than four charged tracks. We also require 
that the ratio of the second to the zeroth Fox- Wolfram moments, R2 = H2/H0 [22], not exceed 0.5. 
R2 is calculated including all detected charged particles and photons. For events with an electron 
in the momentum interval of 2.0 to 2.6 GeV/c, these two criteria reduce the non-55 background 
by a factor of about 6, while they retain more than 80% of the signal events. 

In semileptonic B decays, the neutrino carries sizable energy. In events in which the only 
undetected particle is this neutrino, its energy and direction can be inferred from the missing 
momentum in the event, Pmiss- Wc estimate Pmiss from the difference between the four-momentum 
of the two colliding-beam particles and sum of the four-momenta of all detected particles, charged 
and neutral. To improve the reconstruction of the missing momentum, we impose a number of 
requirements on the charged and neutral particles. Charged tracks are required to have a minimum 
transverse momentum of 0.2 GeV/c and a maximum momentum of 10 GeV/c in the laboratory 
frame. They arc restricted in polar angle to —0.82 < cos 9iab < 0.92 and should originate close 
to the beam-beam interaction point. The detected energy of an individual photon is required to 
exceed 30 MeV. The selection of semileptonic decays can be greatly enhanced by requiring for the 
missing momentum that |pmws| exceed 0.5GeV/c, and that Pmiss point into the detector fiducial 
volume, —0.9 < cos 9miss < 0.9, thereby effectively reducing the impact of particle losses close to 
the beams. Furthermore, since in semileptonic B decays, the neutrino and the electron are emitted 
preferentially in opposite directions, we require that the angle Aa between these two particles 
fulfill the condition cos Aa < 0.4. These criteria on the missing momentum reduce the continuum 
background from QED processes and e+e" qq production by an additional factor of 3, while the 
signal loss is less than 20%. 

The detection efficiencies are estimated by Monte Carlo simulation. With the stated selection 
criteria the efficiency (including effects of brcmsstrahlung) for detecting B X^ev decays is close 
to 50% and largely independent of the electron momentum between 1.1 and 2.0 GeV/c; it gradually 
decreases above 2.0 GeV/c, reaching 35% at 2.6 GeV/c. 

3.2 Background Subtraction 

The spectrum of the highest momentum electron in events selected by the criteria described above 
is shown in Fig. la, separately for data recorded on and below the T(45) resonance. The data 
collected on the T(45') resonance include contributions from BB events and non-BB background. 
The latter is measured using off-resonance data, collected below BB production threshold, and 
using on-resonance data above 2.8 GeV/c, i.e., above the endpoint for electrons from B decays. 
The BB background to the B — >• Xy^ev spectrum is estimated from Monte Carlo simulation with the 
normalization of the individual contributions determined by a fit to the total observed spectrum. 
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1.1 1.5 1.9 2.3 2.7 3.1 3.5 

Electron Momentum (GeV/c) 

Figure 1: Electron momentum spectrum in the T(4S') rest frame: (a) on-resonance data (solid 
circles), scaled off-resonance data (open circles). The solid line shows the result of the fit to the 
non-BB events using both on- and off-resonance data in the interval Pe = 1-1 — 3.5GeV/c; (b) 
on-resonance data after subtraction of the fitted non-BB background (open triangles) compared to 
simulated BB background (histogram) ; (c) on-resonance data after subtraction of all backgrounds 
(open squares), compared to the simulated B — > X^ev events (histogram). The error bars indicate 
statistical errors only. 



3.2.1 Non-BB Background 

To determine the non-BB background we perform a fit to the off-resonance data, collected 

below BB production threshold, and to on-resonance data in the momentum interval of 2.8 GeV/c 
to 3.5 GeV/c. Since the cm. energy for the off-resonance data is lower by about 0.4%, we scale 
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Table 1: Summary of sample composition: principal backgrounds and remaining signal in units of 
10^ events, as well as the signal efficiencies. The errors are statistical, but for the non-BB and 
XcCv backgrounds the uncertainties in the fitted scale factors are included. 



Ap (GeV/c) 


2.0-^2.6 


2.1^2.6 


2.2^2.6 


2.3^2.6 


Total sample 


609.81±0.78 


295.76±0.54 


133.59±0.37 


65.48±0.26 


Non-BB 


142.33±0.63 


105.15±0.48 


74.83±0.36 


50.10±0.25 


XcCV 


423.25±1.58 


160.17±0.86 


39.76±0.35 


4.18±0.09 


J/ip and ip' 


5.50±0.14 


3.55±0.09 


2.07±0.06 


1.04±0.03 


Other e=^ 


1.57±0.04 


0.61±0.02 


0.23±0.01 


0.07±0.01 


TT mis-ID 


1.19±0.04 


0.87±0.03 


0.57±0.02 


0.31±0.02 


K mis-ID 


0.44±0.02 


0.24±0.01 


0.12±0.01 


0.05±0.01 


Other mis-ID 


0.24±0.01 


0.13±0.01 


0.07±0.01 


0.03±0.01 


Other Xj^ei/ 


1.33±0.07 


0.54±0.03 


0.16±0.01 


0.03±0.01 


XuCv signal 


33.97±1.92 


24.50±1.17 


15.80±0.64 


9.67±0.38 


Xucy efficiency (%) 


43.2±0.3 


42.3±0.3 


41.3±0.4 


40.4±0.5 



the lepton momenta by the energy ratio to compensate for the difference. 
The relative normalization for the two data sets is 

TL = 'OFF SLoNdt ^ ^ Q^^^ ^ Q 

SON jLoFpat 

where ,s and L refer to the cm. energy squared and luminosity of the two data sets. The sta- 
tistical uncertainty of is determined by the number of detected pairs used for the jLdt 
measurement; the systematic error of the ratio is estimated to be 0.25%. 
The for the fit to the non-SB events is defined as follows, 

2 _ V- ifid,Pt) -rLUjf ^ ^ ifia,Pj) - Njf 

i 'L^''' j fe>2.8GcV/c) 

Here Nj and nj refer to the number of selected events on- and off-resonance in the j-th or i-th 
momentum bin, and a is the set of free parameters of the fit. f{a,p) is the function approximating 
the momentum spectrum, for which we have chosen an exponential expression of the form 

/(a,_p) = ai +cxp(a2 + a3P + a4P^), (2) 

The fit describes the data well, = 70 for 58 degrees of freedom. Above 2.8 GeV/c, we observe 
(36.7 ±0.2) X 10^ events in the on-resonance data, whereas the fit predicts (36.6 it 0.2) x 10^ events. 

3.2.2 BB Background 

The electron spectrum from S-meson decays is composed of several contributions, dominated by the 
various semileptonic i?-meson decays. Hadronic B decays contribute mostly via hadron misidenti- 
fication and secondary electrons from decays of D, J ftp, and mesons. 

We estimate the total background by fitting the observed inclusive electron spectrum to the sum 
of the signal and individual background contributions. For the individual signal and BB background 
contributions, we use the Monte Carlo-simulated spectra and treat their relative normalization 
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factors as free parameters in the fit. The non-BB background is parameterized by the exponential 
function f{a,pi), as described above. We expand the definition as follows, 



9 _ V- U{a,Pi)-rLnif ^ {f{a,pj) + S{h,pj)-Njf 



where the first sum is for the off-resonance data and the second sum for the on-resonance data. The 
BB electron spectrum is approximated as S{b,pj) = Ylk ^kQkiPj), where the six free parameters bk 
are the correction factors to the default branching fractions for the individual contributions gkipj) 
representing the signal B — > X^ev decays, the background B — > Dev, B D*eu, B — ^ D**ev, 
Trei/ decays, and the sum of other backgrounds with either electrons from secondary 
decays or misidentified hadrons. (Jj mc is a statistical error of the number of simulated events in 
the j-th bin. The momentum spectra gk{Pj) are histograms from Monte Carlo simulations. For the 
largest contribution, the B D*ev decay, we have adjusted the spectrum using the recent BABMi 
measurements of the form factor parameters [20] . 



3.2.3 Fit to Inclusive Spectra 

The fit is performed for the electron-momentum range from 1.1 GeV/c to 3.5 GeV/c, in bins of 
50MeV/c. The lower part of the spectrum determines the relative normalization of the various 
background contributions, allowing for an extrapolation into the endpoint region above 2.0 GeV/c. 
To reduce a potential systematic bias from the assumed shape of the signal spectrum, we combine 
the on-resonance data for the interval from 2.2 to 2.8 GeV/c into a single bin. The lower limit of 
this bin is chosen so as to retain the sensitivity to the steeply falling BB background distributions, 
while containing a large fraction of the signal events in a region where the background is low. As 
the limit is lowered to 2.0 GeV/c, the resulting error on the background subtraction increases. 

The results of the fit and the subtraction of the fitted non-BB and BB backgrounds arc shown 
in Fig. 1, and compared to Monte Carlo simulations. Above 2.3 GeV/c, the non-BB background is 
dominant, while at low momenta the semileptonic BB background dominates. Contributions from 
hadron misidentification are small, varying from 6% to 4% as the electron momentum increases. 
The fit has a of 100 for 75 degrees of freedom. 

Table 1 shows a summary of the data, principal backgrounds and the resulting signal. The 
errors are statistical, but for the non-BB and XcSU background they include the uncertainties of 
the fitted parameters. The data are shown for four different signal regions, ranging in width from 
600 McV/c to 300 MeV/c. We choose 2.6 GcV/c as the upper limit of the signal regions because 
at higher momenta the signal contributions are very small compared to the non-BB background. 
The number of signal events in a given signal interval is taken as the excess of events above the 
fitted background. 



4 Systematic Errors 

The principal systematic errors originate from the fits to the backgrounds, due to the uncertainties 
in the simulated momentum spectra of the various contributions. The uncertainty in the event 
simulation and thereby the impact of the event selection on the momentum dependence of the 
efficiencies for signal and background arc the experimental limitations of the current analysis. In 
addition, there are smaller corrections to the momentum spectra due to variations in the beam 
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energies, and radiative effects. A summary of the systematic errors is given in Table 2 for four 
intervals in the electron momentum. 



Table 2: Summary of the relative systematic errors (%) on the partial branching fraction measure- 
ments for B Xy^ev decays. 
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R — *■ D*]u form f/^rtnr 


1.7 


2.6 


1.7 


0.8 


B Dlv form factor 


0.3 


0.9 


0.8 


0.7 


B — > D**ev spectrum 


2.3 


2.0 


2.6 


1.3 


J/V' and V'' background 


0.8 


0.7 


0.7 


0.5 


Other background 


0.5 


0.2 


0.1 


0.1 


TT mis-ID background 


0.7 


0.7 


0.8 


0.7 


mis-D background 


0.4 


0.3 


0.2 


0.2 


Other hadron mis- ID background 


0.2 


0.2 


0.1 


0.1 


B X^ev background 


1.2 


0.7 


0.3 


0.1 


B momentum 


1.1 


1.5 


1.3 


0.5 


Ngg normalization 


1.1 


1.1 


1.1 


1.1 


Bremsstrahlung and FSR correction 


0.8 


1.1 


1.2 


0.7 


Total Systematic Error 


11.1 


9.8 


8.5 


7.0 



4.1 Detection and Simulation of B — > XuGu Decays 

The detection efficiency for B — > X^^iv decays is determined by Monte Carlo simulation. We include 
in the uncertainty of the signal spectrum not only the uncertainty in simulation of the detector 
response, but also the uncertainty in the simulation of the momentum and angular distributions of 
the charged lepton, as well as the hadrons and neutrinos. 

4.1.1 Detector related uncertainties 

For a specific model of the signal decays there arc three major factors that determine the efficiency: 
the track reconstruction for the electron, the electron identification, and losses due to the detector 
acceptance and the event selection. 

The uncertainty in the tracking efficiency has been studied in detail and is estimated to be 
~ 0.7% per track. The average identification efficiency for electrons above 1.0 GeV/c is estimated 
to be 91.5% [13], based on large samples of radiative Bhabha events and two-photon interactions. 
In BB events the actual efficiencies are slightly lower due to higher track multiplicity. This dif- 
ference decreases gradually from about 2.5% at 1.0 GeV/c to less than 0.8% above 2.0 GeV/c. 
An independent estimate of 0.6% for this uncertainty was derived from a comparison of the effi- 
ciency in data and simulation, for electrons from J/ipK^*^ decays. We adopt an error of 2% for 
the systematic uncertainty in the electron-identification efficiency at 1.0 GeV/c, decreasing to 0.8% 
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above 2.0 GeV/c. In addition, we take into account the impact of the momentum-dependent un- 
certainty of electron identification efficiency on the observed electron spectrum for both signal and 
background (see below). 

4.1.2 Uncertainties in the signal spectrum 

The momentum distributions of the signal electrons are not precisely known because many of the 
B — > Xuii^ decay modes are still unobserved or poorly measured due to small event samples, and 
even for observed ones the decay form factors are not measured. For decays with low mass charmless 
hadrons the simulation is based on ISGW2 model. For decays to higher mass, mostly non-resonant 
states we rely on the model of de Fazio and Ncubcrt [16] and a fragmentation model [17]. 

To evaluate the sensitivity of the signal efficiency to the decay multiplicity and the angular 
and momentum distributions, we randomly vary the individual branching ratios for decays to 
resonant and non-resonant charmless hadrons by 50%, except for B Tr£u and B p£v which are 
currently measured to 30% and 25%, respectively. We observe changes of less than 3.0% for the 
whole spectrum, and less than 1.0% in the signal yield above 2.3GeV/c. 

The systematic uncertainties inherent in the modeling of the signal decays to non-resonant 
hadronic states have been studied by varying the shape-function parameters as determined from 
the measurement of the inclusive photon spectrum [6]. These variations translate to changes of the 
signal lepton spectrum and thus impact the fit to the overall spectrum. The resulting changes in 
the signal branching fraction have been evaluated, taking into account the errors and correlation 
of the shape function parameters. 

Not included in this estimate is the sensitivity to the event selection criteria, specifically those 
based on the variables i?2 and Pmiss- These criteria not only influence the signal, but more so the 
background distributions, and they are discussed below. 

4.2 Non-BB Background 

Systematic errors in the subtraction of the non-BB background could be introduced by the choice 
of the fitting function describing the lepton spectrum and by the uncertainty in the relative nor- 
malization of the on- and off-resonance data. 

To assess the uncertainty in the shape of this background we have compared fits with difi"erent 
parameterizations of the fitting function. In addition to the exponential function described above, 
we have tried linear combinations of Chebyshev polynomials up to the fifth order. The resulting fits 
are equally consistent with the data. The differences in the noii-BB background estimates between 
different parameterizations are at a level of less than 0.5%. 

If relative normalization is treated as a free fit parameter, its deviation from the value based on 
luminosity measurements is less than one standard deviation, which is equal to 1.5%. Thus for the 
normalization the more accurate value based on luminosity measurements is used. As a systematic 
error of non-BB background we take 0.5% of this background contribution which includes the 
errors of normalization factor and background shape approximation. 

4.3 B J/ipX Background 

J/il) and ■0' decays to lepton pairs are vetoed by a restriction on the di-lepton effective mass. 
However, this veto is only about 50% efficient, mostly because of acceptance losses. The remaining, 
mostly single track background is estimated from simulation. We observe a difference of 5.0 it 2.7% 
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between the veto efficiency for lepton pairs in data and simulation, and thus assign a 5% error to the 
residual background. Since this background amounts to 18% (10%) of the signal for Pe > 2.0 GeV/c 
{pe > 2.3 GeV/c) for both differential and integrated spectra, the resulting uncertainty on the signal 
branching fraction is estimated to vary from 0.8% to 0.5%. 

4.4 BB Background 

A major concern in this analysis is the uncertainty in the estimate of the signal events from 
B Xu^v decays obtained from the fit of the sum of various background contributions to the 
observed spectrum. The shapes of BB backgrounds are derived from Monte Carlo simulation. 
The branching fractions for exclusive semileptonic B Xcii' decays are currently not precisely 
known. Thus the lepton spectra from B — > Xcii' decays may differ from those of the simulation. 
For this reason we have introduced scale factors in the fits to the spectrum to adjust the relative 
normalization of the various contributions. To test the sensitivity to the shape of the dominant 
contributions, we have varied the form factor for decays to D*eu and Dev, and changed the relative 
proportion of contributions from narrow and wide resonances to Z)**ez/ decays. 

The differential decay rate for B D*lv can be described by three amplitudes, which depend on 
the three parameters: p^, and i?2- Their measured values are = 0.769 ±0.039 ±0.019 it 0.032, 
Ri = 1.328 ± 0.055 ± 0.025 ± 0.025, and i?2 = 0.920 ± 0.044 ± 0.020 ± 0.013 [20]. The quoted 
errors are statistical uncertainties from the data and Monte Carlo samples, and systematic errors, 
respectively. The B D£u differential decay rate can be described by a single parameter p^. To 
study the impact of form factor variations we reweight the Monte Carlo-simulated spectrum for a 
given decay mode with the relative change of the generator-level spectrum due to changes in the 
form-factor parameters, and repeat the standard fit to the data. From the observed changes in the 
signal rate as a function of the choice of the form-factor parameters for 0*61/ decays, we assess the 
systematic error on the signal rate by taking into account the measured form-factor parameters, 
their errors, and their covariance matrix [20]. For Deu decays, we rely on a measurement by the 
CLEO Collaboration [23], Fd{w)/Fd(1) = 1 - pI>{w - 1), where p% = 0.76 ± 0.16 ± 0.08 is the 
linear slope of the u)-dependence. The variable w is the product of the four-vector velocities of the 
B and D mesons and corresponds to the relativistic boost of the D meson in the B rest frame. 
To estimate the impact of the uncertainty in pj^ we likewise compare the default fit results with 
fits performed with a rcwcightcd lepton spectrum. In this case, we replace the default simulation 
based on the ISGW2 model with the simpler form factor formulation and vary the measured value 
of pI, by one standard deviation. We adopt the mean shift of the signal rate as a systematic error. 

To assess the impact of the poorly known branching fraction for various D**ev decay modes 
on the shape of the lepton spectrum we have repeated the fit with the different relative branching 
fractions for the individual decays modes. As long as we do not eliminate the decays to the 
two narrow resonances, £)i(2437) and £)2(2459), we obtain reasonable results. Specifically, if we 
eliminate the decays involving the two wider resonances, -Do(2308) and -D^(2460), the results change 
by less than 3%. We adopt this change as the estimate of the systematic error due to the uncertainty 
of decays to D** states. 

Similarly, we vary the branching ratios for secondary electrons from semileptonic D decays by 
10% and adopt the observed change as a systematic error. There is a small background from events 
which contain a B ^ Xuev decay but contribute to the background rather than the signal. We 
estimate the uncertainty in this contribution to be 30%. 

For background from hadronic B decays, the uncertainty in the spectrum is primarily due to the 
uncertainty in the momentum-dependent hadron misidentification. The uncertainties of misidenti- 
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fication probabilities are estimated to be 20%, 30%, and 50% for pions, kaons, and protons respec- 
tively. The uncertainty in the fractions of pions, kaons, and proton is taken as a difference between 
simulated and observed charged particle spectra, which is about 5% for pions and kaons, and 50% 
for much smaller contribution from protons and antiprotons. With these uncertainties in the hadron 
misidentification backgrounds, the fractional error in the number of subtracted background events 
is ~ 20% for pions, ~ 30% for kaons, and 70% for protons. 

4.5 Uncertainty in the B Meson Momentum Spectrum 

The non-zero momentum of the B mesons in the T{4S) rest frame affects the shape of the elec- 
tron spectra near the endpoint. To estimate the systematic error of the inclusive lepton spectra 
associated with the uncertainty of initial B meson momentum spectrum we compare the simulated 
and measured energy spectra for fully reconstructed charged B mesons. The widths of the energy 
distributions agree well, but in some of the data sets we observed a shift in the central value of up 
to 2.2 MeV relative to the simulation. We correct the simulation for these shifts, and assess the 
effect of the uncertainty of 0.13 MeV in this shift. 

4.6 Bremsstrahlung and Radiative Corrections 

For comparison with other experiments and with theoretical calculations the signal spectrum re- 
sulting from the fit is corrected for bremsstrahlung in the detector and for final-state radiation. 
Corrections for QED radiation in the decay process are simulated using PHOTOS [24]. This 
simulation includes multiple-photon emission from the electron, but docs not include electrowcak 
corrections for quarks. The accuracy of this simulation has been compared to analytical calcula- 
tions performed to 0{a) [24]. Based on this comparison we assign an uncertainty of 10 - 15% to 
the PHOTOS correction, leading to an uncertainty in the signal yield of about 1%. 

The uncertainty in the energy loss due to bremsstrahlung is determined by the uncertainty in 
the thickness of the detector material, estimated to be (0.0450 zt 0.0014)Xo at normal incidence. 
This thickness was verified using Bhabha scattering as a function of the polar angle relative to the 
beam. The impact of this uncertainty on the signal rate was estimated by calculating the impact 
of an additional 0.0014Xo of material; it is of the order of 1%. 

4.7 Sensitivity to the Event Selection 

We have checked the sensitivity of the fits to the lepton spectrum to changes in the event selection, 
as well as the momentum dependence of the electron selection efHciencies. These variations of the 
cuts change the signal efficiency and lead to large variations up to 50% in the size of the non-BB 
background and up to 20% in the BB background. Though some of the observed changes in signal 
yield may already be covered by the form factor and other variations, we conclude that these tests 
do reveal significant changes that have to be accounted for. We interpret the observed changes as 
inadequacies in the simulation of the Monte Carlo-simulated lepton spectra and adopt the observed 
changes between the default fits and the looser selection criteria as systematic errors. 

The largest variation (5%) is observed for changes in the restriction of the Fox- Wolfram ratio 
R2 from the default value of 0.5 to 0.4 and 0.6. Other sizable variations are observed for changes in 
the absolute value and direction of the missing momentum vector. R2 and the missing momentum 
are quantities that are derived from the measured momenta of all charged and neutral particles 
in the event and are therefore sensitive to even small differences in data and the simulation. We 
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consider the chosen variations of the cut variables representative for the error estimation, add the 
observed changes in quadrature, and include them in the overall systematic error. 



5 Results 

5.1 Determination of Partial B X^eu Branching Fraction 

For a given interval Ap in the electron momentum, we calculate the inclusive partial branching 
fraction B — > XuCv according to 

Here Ntot refers to the total number of electron candidates detected in the on-resonancc data and 
Nf)g refers to the total background, from non-i?i? and BB events, as determined from the fit to 
the spectrum. e{Ap) is the total efficiency for detecting a signal electron from B Xucu decays 
(including bremsstrahlung in detector material), and 5rad accounts for the distortion of the electron 
spectrum due to final-state radiation. This is a momentum-dependent correction, derived from the 
Monte Carlo simulation based on PHOTOS [24]. The total number of produced BB events is 
iV^^ = (88.36 ± 0.02 ± 0.97) x 10^ 

The differential branching fraction as a function of the electron momentum in the T{4S) rest 
frame is shown in Fig. 2, fully corrected for efficiencies and radiative effects. Partial branching 
fractions for the four different momentum intervals arc summarized in Table 3. There is excellent 
agreement with the preliminary BABAR measurement [11] and also with a measurement by the 
CLEO Collaboration [10]. 



Table 3: Preliminary results on the partial {AB) and the total branching fraction {B) for inclusive 
B X^eu decays for four momentum intervals. /„(Ap) is a fraction of the lepton spectrum, that 
falls into the Ap momentum interval. The values of fu{Ap) were derived by the CLEO Collaboration 
[10] from the shape function parameters based on the measurement of the b ^ sj spectrum. 



Ap(GeV/c) Ag(10-^) /„(Ap) g(10-^) IKhKlO"") 

2.0^2.6 0.480 ±0.029 ± 0.053 0.278 ± 0.043 ± 0.030 1.73 ± 0.22 ± 0.33 3.94 ± 0.25 ± 0.37 ± 0.19 

2.1 H- 2.6 0.355 ±0.018 ±0.035 0.207 ± 0.037 ± 0.027 1.71 ± 0.19 ± 0.38 3.93 ± 0.22 ± 0.43 ± 0.19 

2.2^2.6 0.235 ±0.010 ±0.020 0.137 ±0.025 ±0.023 1.72 ± 0.16 ± 0.42 3.93 ± 0.19 ± 0.48 ± 0.19 

2.3^2.6 0.149 ±0.006 ±0.010 0.078 ± 0.015 ± 0.016 1.90 ± 0.15 ± 0.52 4.14 ± 0.16 ± 0.56 ± 0.20 



5.2 Extraction of the Total Charmless Branching Fraction and \Vub\ 

To determine the charmless semileptonic branching fraction B{B X^eu) from the partial branch- 
ing fraction AB{Ap), one needs to know the fraction fu(Ap) of the spectrum that falls into the 
momentum interval Ap. The CLEO Collaboration [10] has used the measurement of the inclusive 
photon spectrum from b ^ s'j transitions to derive the parameters describing the shape function 
to determine the fraction /„(Ap), and to estimate the uncertainties of the resulting shape function 
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Figure 2: The differential branching fraction for charmless semileptonic B decays (data points) as 
a function of the electron momentum (in the T(4S') rest frame) after background subtraction and 
corrections for bremsstrahlung and final state radiation, compared to the Monte Carlo simulation 
(histogram). The errors indicate the statistical errors on the background subtraction, including the 
uncertainties of the fit parameters. 



parameters. The /^(Ap) fractions arc listed in Table 3, together with the total branching frac- 
tions. Contrary to the CLEO analysis, we have fully corrected the AB for radiative effects, and 
consequently the fractions listed here do not include such corrections [25]. The quoted errors 
on B and /„ arc statistical and systematic. The systematic error of /„ includes the systematic 
error of BB subtraction, uncertainties of the shape function approximation, dependence on as 
scale, and the uncertainty of theoretical prediction of from b ^ sj shape function. The errors 
on B and \Vui,\ are the total error from AB measurement and the error from J^. The last error 
in \Vub\ is the theoretical uncertainty of the translation from B to \Vub\- Surprisingly, the overall 
precision does not depend very strongly on the chosen momentum interval. While the experimental 
errors are smallest for the interval from 2.3 to 2.6 GeV/c, the dominant uncertainty arises from the 
determination of the fraction and this increases substantially with higher momentum cut-off. 
Thus we quote as the branching fraction measurement the result based on the data in the interval 
2.0-2.6 GeV/c, 

B{B ^ Xueu) = (1.73 ± 0.22e^p ± 0.33/J x lO'^, (5) 

where the first error represents the sum of the statistical and systematic error on the partial 
branching fraction, and the second error refers to the uncertainty in the measured photon spectrum 
and the extraction of the shape function. Prom the inclusive charmless semileptonic branching 
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fraction and the average B lifetime we can extract V^b) 

1 /2 

= 0.00424 (^^i^^l^l:^^ (1.0 ± 0.028oPi. ± 0.039„J. (6) 

Here, we rely on a formulation of [26] and [27], taking into account recent measurements of the 
B lifetime of 1.604 ± 0.012 ps [28], the b quark mass in the kinetic mass scheme, m(,(lGeV) = 
4.61 ± 0.07 GeV, and other parameters of the heavy quark expansions [29]. The first error represents 
the linear sum of the perturbative and non-perturbative QCD corrections, the second error is due 
to the uncertainty in m^. An overall correction of 0.7% is included to account for QED corrections. 
The results for the four momentum intervals are presented in Table 3. They are consistent with 
each other. For the interval from 2.0 to 2.6 GeV/c we quote as a preliminary result, 

IK&I = (3.94 ± 0.25e^p ± 0.37/„ ± 0.19*^2,) x lO'^. (7) 

Here the first error represents the total experimental uncertainty, and the second refers to the 
uncertainty on the determination of the fraction /„ from the B Xgj decays (taken from the 
CLEO analysis), and the third combines the stated theoretical uncertainties in the extraction of 
|V^b| from the branching ratio. No additional uncertainty due to the theoretical assumption of 
quark-hadron duality has been assigned. 

Recently new preliminary determination of the shape function parameters based on analysis 
of 6 ^ S7 spectrum by the Belle collaboration became available [30] with the central values of 
a = 2.27 and m^^ = 4.62 GeV/c^ . The Belle shape-function parameters differ from those of CLEO 
by more than one standard deviation. The results for B — > Xucv partial and total branching ratios 
based on this new set of shape-function parameters and their errors are presented in the Table 4. 
The quoted errors on AB are statistical and systematic, respectively. The total error on /„ includes 
the statistical and systematic errors of the shape-function measurement and the uncertainty of the 
theoretical calculation of fu from b ^ sj shape function, estimated as in [10]. The errors on B and 
\Vub\ are the total error from AB measurement and the error from /„, respectively. The last error 
in \Vub\ is the theoretical uncertainty of the translation from B to [V^fej. 



Table 4: Preliminary results on the partial (AB) and the total branching fraction {B) for inclusive 
B Xy^ev decays for four momentum intervals. The recent Belle measurement of 6 ^ 57 spectrum 
was used [30]. 

Ap(GcV/c) Ag(10-^) /„(Ap) g(10-^) \Vub\{lO-'') 

2.0 -=-2.6 0.531 ±0.032 ±0.049 0.246 ±0.031 2.16 ± 0.24 ± 0.27 4.40 ± 0.24 ± 0.28 ± 0.21 
2.1 2.6 0.381 ±0.020 ±0.034 0.174 ±0.026 2.19 ± 0.23 ± 0.33 4.44 ± 0.23 ± 0.33 ± 0.21 
2.2 -=-2.6 0.245 ±0.011 ±0.020 0.110 ±0.022 2.23 ± 0.21 ± 0.45 4.47 ± 0.21 ± 0.45 ± 0.22 
2.3 H- 2.6 0.153 ±0.006 ±0.011 0.058 ±0.017 2.64 ± 0.22 ± 0.77 4.87 ± 0.20 ± 0.71 ± 0.23 



In conclusion, we have a preliminary measurement of the differential lepton spectrum for charm- 
less semileptonic B decays above 2.0 GeV/c, and have extracted the CKM matrix element \Vuh\ with 
improved experimental accuracy. Further improvements to this measurement are expected from 
new measurements of the inclusive photon spectrum from B — * Xg^ decays and further advances 
in our theoretical understanding of the shape functions, their relation to the parton-level lepton 
spectrum, and effects of the fragmentation of the s and u quarks. 
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